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Abstract 

A survey was made of all interplanetary shocks detected by the plasma 
science experiment aboard the Voyager 1 spacecraft between 1.2 and 8.5 AU. 
Shock nox-mals and shock velocities are determined. The variation of shock 
frequency and various shock parameters with heliocentric distance is 
discunsed. The following results are observed: l) Beyond 1.2 AU, the vast 

majority of shocks were associated vith interaction regions between high and 

»’ • * * , 

low speed streams; of 95 events, only 1 was clearly associated with a 
transient event. 2) Forward shocks were more numerous and seemed to form 
closer to the sun than reverse shocks. 3 ) Forward shocks were stronger than 
reverse shocks. 

The energy balance of three shocks is examined. A close agreement is 
found between the measured and the predicted pressure ratios across these 
shocks. The contribution of shocks to the global energy balance is discussed. 
Shocks are found to have a significant effect in heating the solar wind. 
However the effect of heating due to shocks appears to be inadequate to 
explain the heat flux in the solar wind between 1 and 9 AU. 
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Introduction 

An understanding of the origin, evolution, and eventual fate of Bhocks in 
the outer heliosphere is important if one is to understand the dynamics and 

t 

energetics of the solar wind beyond 1 AU. 

Observations of interplanetary shocks near 1 AU have been reported by 
numerous observers over the past two decades [Sonnet et al, 1964; Dryer et al, 
1972, 1975? Bavassano et al, 1973} Hundhausen et al, 1970; Chao and Lapping, 
1974]* More recently observations at larger heliocentric distances have been 
reported using data from Pioneer 10 and 11 [Smith and Wolfe: 1976, Smith and 
Wolfe; 1977; Mihalov and Wolfe, 1979]. 

Shocks in the solar wind can originate from flares and other solar 
transients, or they can be caused by the interaction between high and low 
speed streams of the solar wind. Moat shocks in the inner heliosphere (R < 1 
AU) can be associated with transient events [Chao and Lepping, 1974]» 

Extensive theoretical and experimental work has been done on the modeling of 
flare-generated shocks [Hundhausen and Gentry, 1969; Stoinolfsen et al, 1975; 
Borrini et al, 1982]. 

Beyond 1 AU, the principle source of shocks appears to be stream 
interactions [Smith and Wolfe, 1979]* This effect was predicted by Sarabai 
[1963], Parker [1963], Colburn and Sonnet [1966], and others. The name 
'interaction region' was coined for the plasma affected by these stream 
interactions by Smith and Wolfe [1976], These regions have since come to be 
called 'co-rota ting interaction regions' (CIR's) since the stream structure 
persists for several rotation periods and the entire quasi-stationary 
structure rotates with the solar rotation period. The configuration of low 
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speed stream, forward shock, shocked plasma, reverse shock, and high speed 
stream is diagrammed in Figure 1. 

Numerical models for the formation and evolution of CIR's have been 

* 

constructed UBlng a two-dimensional MHD approximation [Hundhausen, 1973a, 
1973b; Dry^r et al, 1978; Steinolfsen and Dryer, 19781. A comparison between 
theory and observation has beeen made for the limited periods during multiple 
spacecraft line-ups [Burlaga et al, 19831. More recently, models have been 
devised. using, a full set of three-dimensional MHD equations {Pizza, 1983). 
However no attempt has been made to compare model predictions concerning shock 
evolution with observations of a large (>12) sample of shocks. Furthurmore, 
none of these models treat shocks explicitly as such. Instead, finite 
dissipative terms are introduced the handle regions containing a large Jump in 
plasma parameters. 

In order to draw reliable conclusions concerning the evolution of shocks 
in the outer heliosphere, one must conduct a survey of a large number of 
shocks. To learn the average contribution of shocks to heating of the solar 
wind, one must have a sufficiently large quantity of shock data that the 
averages are statistically meaningful. Finally, a new study of interplanetary 
shocks between 1 and 10 AU using Voyager spacecraft data should provide a 
useful comparison with Pioneer observations in this regions because the 
Voyager and Pioneer spacecraft traversed this region in very different 

N v. 

portions of the solar cycle. ‘ 

The Voyager 1 spacecraft was launched in late 1977* A plot of the 
Voyager spacecraft trajectories together with the orbits of Earth, Jupiter and 
Saturn is shown in Figure 2. Between day 300 of 1977 and day 262 of 1980, the 
spacecraft traversed the region between 1.2 and 8.9 AU from the 3un. During 
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this time, a large number of shock events in the solar wind were observed by 
the plasma science experiment aboard the spacecraft. Due to the irregular 
nature of the solar wind stream structure during the rising portion of the 
solar cycle when these observations were made, it is sometimes difficult to 

i 

determine which features are transient events and which are co-rotating, 
especially if flares arise from longitudes near to the source regions of high 
speed streams. Nevertheless, only one of the shocks that were observed was 
associated with a stream feature that did not recur for at least three solar 
rotations: the shock of day 260 of 19T9* 

Background 

Co-rotating shocks are formed by the interaction of high and low speed 
Btreams in the solar wind. In principle, four shocks could form as a result 
of the interaction. Two would be fast and slow forward shocks, propagating 
downstream into the low speed stream, and two would be fast and slow 'reverse' 
shocks, which are convected downstream while propagating upstream into the 
high speed stream. Usually the two slow shocks are not observed. Between the 
forward and reverse shocks a region of shocked gas of intermediate velocity 
and enhanced density and temperature forms. In the Voyager 1 data, 
stream- produced shocks appear only beyond a heliocentric distance of 2-3 AU. 
Thereafter stream-produced shocks appear frequently with one or more shocks 
appearing per solar rotation. 

At a shock, the MHD equations become the Rankine-Hugoniot relations 
describing the Jumps in the plasma and the magnetic field parameters. If the 
bracket operator, I), represents the difference between the post- and 
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pre-shock values , the subscripts n and t represent the components normal to 
and parallel to the shock, and U is bulk velocity in the shock rest frame, the 
relations become: 


lpu n l = 0 

mass flux 

(1) 

= 0 

tangential Bhear balance 

(2) 

[pU^+p+Bf/STt) = 0 
n u 

normal pressure balance 

(3) 

(B n ) = 0 

normal B-field continuity 

(4) 


tangential E-field continuity 

(5) 


where p is the mass density (gm cm” 2 ), and p is the isotropic pressure. 

Assuming a polytrope approximation, one can write an energy conservation 
equation: 


|i(uW) + -L £ + ri-B? ♦ -j^l = 0 

2 t n y- 1 p 4itp t 4npU 


. . 0 Bo. B U 

1 *>2 t n n, 


( 6 ) 


However there is no reason to assume that the polytrope approximation is 
valid in the region of a shock. 

From equations l) , 2), and 4), one can obtain: 


lB t (f mVV U * )] = 0 


where the normal mass flux f = pU_ 

m r n 


(7) 


T 


One can then define three possible cases: 

Case 1: f U, -B^/hu > 0 Fast Shock 

m In n 

In the case of a fast shock, the Alfven mach number of the shock is 
greater than one. Since f and are continuous, the Jump in density across 
the Bhock has the same sign as the Jump in the B~ field magnitude. 

Case 2: fJU., -B*jArc < 0 Slow Shock 

m in n 

In the case of a slow shock, the Alfven mach number is less than one, and 
the Jump in density across the shock has the opposite sign from the Jump in 
B- field magnitude. Since the shock speed is less than the Alfven speed, a 
alow shock could in principle decay into a pressure wave. 

O 

Case 3a: f„U, -B^A* = 0 Rotational Discontinuity 

m in n 

A, rotational discontinuity, or ’Alfven shock’ is essentially a finite 
amplitude Alfven wave. The plasma bulk velocity does not vary across the 
discontinuity and the discontinuity propagates at the Alfven speed relative to 
the unshocked plasma. 

Case 3b: U n = 0 Tangential Discontinuity 

Tangential discontinuities do not propagate in the plasma. Tangential 
discontinuities frequently occur at stream interfaces [Gosling et al. 1978). 
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Unlike shocks, neither rotational nor tangential discontinuities 

* 

transform bulk kinetic energy into thermal energy. Thus , neither fora of 
discontinuity affects the energy balance. Neither discontinuities nor their 
consequences are studied in this paper. 

The Method of Locating Shocks 

An attempt was made to locate all of the shocks observed T>y the Voyager 1 
spacecraft. No similar attempt has yet been made to study the Voyager 2 data 
due to the fact that the Voyager 2 data are noisier and have many data gaps. 
Several competing methods were used to search the Voyager 1 plasma data for 
shocks. 

In the initial search, 27-<J*y plots were made of hourly averages of the 
Voyager 1 plaBma data from Day 26 q 1977 to fay 300 I98O. The density and 
temperature profiles were visually examined for correlated Jumps in density 
and temperature. Approximately 70-80^ of the shocks observed were detected by 
this visual search. The remaining 20-30^ of the shocks were too small to be 
stand out prominently. 

In the second search, a computer algorithm was written to list all 
velocity Jumps of greater than 20 km-sec - ^ in the Voyager 1 hourly averages 
which had correlated or anti-correlated density ’Jumps. Approximately 98^ of 
the shocks were detected by the computer search, but it picked up 2-3 spurious 
events for every real shock. These spurious events were due to data gaps , 
slow variations of solar wind parameters which appeared discontinuous in the 
hourly averaged data, or fluctuations in velocity which proved not to be 
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correlated with Jumps in density and temperature when the high resolution data 
were examined. These causes were all of roughly equal importance. 

These two lists were then compared with plots of the highest resolution 
plasma data from Voyager 1. Table 1 shows a final list of those events which 
might be candidates for shocks. Events which involved a discontinuous Jump in 
the plasma data were listed as 'shocks'. Events which involved changes in the 
plasma data which resembled shocks but had a finite width of several minutes 
were listed as 'slopes'. Data gaps across which changes in plasma bulk 
velocity, density, and temperature occured which were consistent with the 
occurrence of a shock were called 'gaps' . The lengths of these gaps are 
indicated in Table 1. 

An independant search was made of three months of the plots of the high 
resolution plasma data between 3 and 4 AU to determine how many shocks might 
have been missed by the search of the hourly average data. No new shocks were 
located . 

Due to the presence of waves and random fluctuations, it is difficult to 
search for shocks with small velocity Jumps. Slow variations of solar wind 
parameters frequently cause small apparently discontinuous Jumps in the hourly 
averaged data. For this reason, the hourly averaged data were only searched 
for velocity Jumps of greater than 20 km-sec** 1 . Many of the shocks located by 
these searches proved to have velocity Jumps of less than 20 km-sec" 3- when the 
high resolution data were examined. But these surveys still have a lower 
detection threshold of approximately 20 km-sec" 1 in velocity Jump. 

Figure 3 shows a histogram of the shock velocity jumps. This histogram 
is strongly peaked near the detection threshold. It is unclear whether this 
peak is due to a real absence of shocks with velocity Jumps of less than 20 
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km-sec - ^, op whether It is due to selection effects. However a more recent 
search of a limited period of Voyager 2 spacecraft data suggests that there 
are a large number of shocks lying below this detection limit. 

The Method </: 4 ' Determining Shock Normals 

« 

A large number of techniques have been proposed for determining shock 

normals [Colburn and. Sonnet, 19 66 ; Lepping and Argentiaro, 1971 1 . The method 

• • < • , . ’ *■ • 

used here utilizes both plasma and magnetic data, but in the interest of speed 
it does not attempt a least squares fit to all of the Rankine-Hugcmiot 
relations. 

From Equations 4 and 5 , one can obtain the r coplanarity theorem* [ Colburn 
and Sonett, 1966] which states that the shock normal must be aligned with: 

( 8 ) 


where the subscripts 1 and 2 refer to the values of the unshocked and the 
shocked parameters respectively. 

There are several practical problems in the use of this expression: The 

pre- and post-shock magnetic field are often nearly parallel, particularly for 
qusai-perpendicular shocks ( shocks with normals approximately perpendicular to 
the B-field in the unshocked plasma) . Because of fluctuations in the post 
shock magnetic field, it is sometimes difficult to determine with precision. 
These problems can be mitigated by utilizing the velocity data: From the mass 

flux equation and the continuity of the normal B-field one can obtain an 
alternate form of the co-planarity theorem [Abraham-Sehrauner and Yun, 1976] : 
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n II where A$ a (Jg-fc^/ltg-S^ io a unit vector (9) 

The uncertainty in this expression is generally much smaller than that in 
Equation 8 except for the case of quasi-parallel shocks. 

If the magnetic field is small so that the dynamic pressure is very much 
greater than the dynamic pressure (M^>>1, such as is the case for faBt 

shocks), then a third expression for the shock normal can be derived from 

* , 

Equation 1. One can obtain the following expression for the shock velocity 
(Lepping and Argent iero, 1971 h 

* s = [i'(p 2 V 2 -P 1 V 1 )/(p 2 -p 1 )]n (10) 

If one neglects the magnetic pressure in the Rankine-Hugoniot relations, the 
tangential shear balance equation becomes: 

IpU n u t ) so (U) 

But (ll) and (l) together imply that (V^.1 2 0 across the shock, or: 

(p a Wi> <„, 

s p 2 - Pl 

Equation 12 is rigorously true if B=0. 

The error due to neglecting the contribution of the magnetic shear to the 
tangential shear balance equation is of the order of the ratio of the two 
terns in Equation 2, or: 



For fast shocks, the shock normal is calculated using (12) and compared 
with the normal obtained from the alternate co-planarity relation (9) . If the 
two normals agree, the normal determination is considered to be accurate. For 
slow shocks the shock normals and velocities are determined using only 
Equation 9. 

* , » ' * 

A typical event, the forward shock seen at 0455 day 313 1979* is shown in 
Figure 4. The various shock parameters are measured for the data point 
immediately preceeding and immediately following the shock. If significant 
oscilations exist in some parameter, the peak values of the oscillation 
closest to the shock are measured. The mean value and half of the 
peak-to-peak value of that oscillation are taken to be the mean value and the 
uncertainty respectively of that parameter. Other methods for the estimation 
of shock parameters are possible, this particular method was chosen for 
reasons of convenience. 

Typical fluctuations in the radial velocity are on the order of 1$. The 
fluctuations in the other parameters are larger, usually of order 10$ and 

t 

sometimes as large as 30$. For quasiperpendicular shoci's, this implies an 

angular uncertainty in the shock normal determination of approximately 15° and 

-1 

an uncertainty of at least 10 km-sec in the velocity determination. 

Quality of the Determinations of Shock Normals 
Two techniques are used to evaluate the quality of the shock normal 
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determination. The first, mentioned in the proceeding section, is to compare 

* 

the direction of the shock normal determined through use of the continuity 

equation with the direction of the shock normal determined using the 

coplanarity theorem. This technique is applicable only for fast shocks* 

Figure 5 shone a histogram of the angle between the normals determined 

for each fast shock by the two method?. The largest source of disagreement 

is due to the error in neglecting the effect of the magnetic pressure when 

UBing cnly the continuity equation to determine the shock normal. This error 

is expressed in Equation 13. If is of the order 3, then 60 is typically 
o 

of order ±20 . For 64 of the 71 fast shocks observed, the disagreement 
between the two normal determinations was of this order. The remaining 7 
cases were examined in detail. In each case, the disagreement was 
attributable to unusually large uncertainties in the measurements of the basic 
shoqk parameters. ' 

tf'rom the mass flux equation, the tangential E-field continuity equation, 
and the tangential shear balance equation, one can obtain: 


N 1 3l t|W B 2t , w-2i 

M- 0 1+1 B 1 ^ M AnJ 

n 2 a 2t a it An 


where the perpendicular Alfven Mach number M 


‘An 




The uncertainty in (l4) can be written: 


(l4) 


N 6B.. B„. 6B-. N. B 0 . _ 6M.„ B of _ 

+ #£» - 2 nM 


(15) 


If there are fluctuations of lOj in the measured densities and B-fields, 
it is possible for the predicted and measured density ratios to disagree by a 
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factor of 1.7* For a Bhock with a large Alivenic Mach number, 

M^>('Bg t /B^-l) , the uncertainty in the perpendicular Mach number can be an 
even larger source of error: As the angle 0 between the shock normal and the 

B-field approaches 90°, small uncertainties in 0 can lead to large 
uncertainties in the predicted density ratio. 

Figure 6 shows a histogram of the theoretical density ratio divided by 
the density ratio measured for each shock. The normal direction was 
determined using the form of the co-planarity relation in Equation 9. For 69 
of the 82 events, the agreement of the two density ratios is within a factor 
of two. As shown below, this is as good as can be expected given the 
uncertainties in the measured data. 

For all of the events observed at Voyager 1, the error in the 
determination of the shock normal inferred from the error In the predicted 
density ratio is within the uncertainty in the shock normal set by the 
uncertainty in the measurement of the data. For 69 of these 82 events, the 
determinations of the shock normals are accurate to within at least 15°. 

Shocks in the Solar Wind: Observations 

Figures 7 through 10 [Gazls and Lazarus, 1983 1 show the variation of 
various shock parameters with heliocentric distance. 

Figure 7 shows a plot of the frequency of occurence of forward shocks and 
reverse shocks versus heliocentric distance. The frequency of occurence of 
shocks was determined by taking a running average of the time spacing between 
shocks ror groups of twelve shocks with a slip of one shock between successive 
averages. The first two shocks In Table 1 were not included in these plots 
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since they are imediately followed by the Voyager 1 data gap between day 320 
1977 and day 100 1978. 

As can be seen, there is a distinct difference between the radial 
profiles of the frequency of occurence of forward shocks and reverse shocks . 

The forward shocks originate earlier in the solar wind, occur more frequently, 

* 

and continue to be observed out to greater distances than the reverse shocks . 
The fact that the observed shock frequency decreases with increasing 
heliocentric distance does not necessarily mean that some of the shocks 
actually disappear as they propagate outwards from the sun. The shocks might 
continue to exist but drop below our detection threshold of AV=20 km-sec _1 . 

It also probable that the apparent radial variation of shock frequency is in 
fact a temporal variation. In view of the lack of any obvious change in the 
stream structure seen by IMP 8 at 1 AU, and the lack of any change in the 
average of the bulk parameters seen at 1 All, this source Beems unlikely, but 
until a complete list is made of all shocks observed by the Voyager 2 
spacecraft, it will be difficult to separate radial and temporal variations. 

Superimposed upon the radial variation of shock frequency are large 
temporal variations, particularly for the forward shocks in late 1978 when the 
Voyager 1 spacecraft was near 5.5 AU. These temporal variations do not 
correspond to any obvious variation in the average plasma bulk parameters 
observed at the IMP 8 spacecraft during the corresponding time, (it is 
interesting to note, though, that the large drop in frequency of forward 
shocks observed at 5.5 AU occurs in mid-late 1979, at the onset of a sudden 
decrease in geomagnetic storm activity) . 

It is instructive to plot the radial variation of various parameters 
related to the strength of shocks; 
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Figure 8 shows a plot of the running average of the speeds of twelve 
successive shocks (forward shocks and reverse shocks) plotted versus 
heliocentric distance. Twelve-shock averages were chosen as the test 
compromise between poor time resolution and excessive fluctuations in the 
average. As' with shock frequency, the shock speeds for forward shocks are 
uniformly higher than the Bhock speeds for reverse shocks . 

Figure 9 shows a plot of a twelve-shock running average of the shock 
velocity Jump for forward shocks and reverse shocks. ‘Since the shock velocity 
jump is related to shock speed, it is no surprise that the velocity Jumps 
across forward shocks are higher and change more slowly with heliocentric 
distance than do the velocity Jumps across reverse shocks. The average of the 
velocity Jump across the reverse shocks declines approximately twice as fast 
as the average of the velocity Jump across the forward shocks. 

Figure 10 shows a plot of twelve-shock running averages of the density 
ratios across the shock for forward shocks and reverse shocks versus 
heliocentric distance. TJnlike the shock frequency, speed, and velocity Jump, 
the density ratio across both the forward and the reverse shocks remains 
fairly constant with increasing heliocentric distance. The twelve-shock 
average of the density ratios across the forward shocks is between 1.6 and 
2.8, though individual shocks have density ratios which range from £1 to 4. 

The reverse shocks are somewhat weaker with an average density ratio of 
approximately 2, but the range of density ratios for individual reverse shocks 
is much the same as for forward shocks. 
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A Model of Shock Evolution 

* 

For a perpendicular fast Bhock with a radial normal, the shock speed is 

* 

determined only by the continuity equation. The equation for shock speed, V a , 
in the spacecraft frame is (12). For purely radial velocities, this can be 
rewritten: 


V s = 


*V V 1 

x-1 


ro 

where the density ratio x = — 

pi 


( 16 ) 


The geometry of forward and reverse shocks as Been in the spacecraft 

i 

frame is shown in Figure 11. For the forward shock, P 2 > P^ and V g >V 2 >V^. For 
the reverse shock, p 2 >p^ and V.j>V 2 >V 8 . 

Taking the convective derivative, (l6) becomes: 


£_y 

Dts 


= 

X-1 


V 2 (x-l)-(xV 2 -V 1 ) D 


X-1 


Dt 


x + X - T=V 


D_ 

Dt' 


i) 


(17) 


The density ratio, x, and the post shock velocity, Vg, vary only slowly 
with respect to R and can be treated as constant as the shocks propagate 
outwards from the sun. Therefore, for typical values of shock parameters such 
as in Table 2, the first two terms on the right side of (IT ) are negligable 
and one can write: 


ILj/ 

Dts 


— i— JLy 

x-1 dTi 


(18) 


The pre-shock velocity V 1 changes due to the fact that the shock Is 
propagating into a velocity gradient. If one takes the convective derivative 
of in the shock frame, one obtains: 



(But the convective derivative in the spacecraft frame vanishes, or:) 


ifo ■ fe*i + v i ki • 0 


( 20 ) 


Using (20) to eliminate the partial time derivative, (19) becomes: 


Dt V l = ^ V s" V l^R V l = U 10R V 1 


( 21 ) 


And (l8) becomes: 


u — v 

Fb x-i 1 airi 


( 22 ) 


or: 


Dll 



U 


1 5H 1 


(23) 


If the shock decays when the upstream flow velocity in the shock frame 
drops below the AlfVen speed, then the lifetime of a fast shock will be 
approximately: 



(24) 


Alternatively since U^fx/x-ljAV, if the shock drops below the detection 
threshold when the velocity jump, AV, drops below some value, AV^., then the 
shock will have an apparent lifetime: 
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\ - ( u i - (25) 

So in the case of a perpendicular fast shock, a shock which initially has 
a larger sh6ck speed will survive longer than a shock which starts with a low 
shock speed. In the case of shocks having the same shock speed a forward 
shock, since it is convected outwards from the sun at a higher velocity than a 
reverse shock, will propagate out to a larger heliocentric distance before it 
disappears . 

Typical numbers for interplanetary shocks at 1+AU are shown in Table 2. 
There is reasonable agreement with the radial profiles observed in Figure 9» 

A Model of Shock Formation 

If one takes the values of the plasma parameters before and after the 
interaction region and attempts to use the six Rankine-Hugoniot relations, 
(l)-(6), to determine -a set of shock velocities and the density ratios which 
will reproduce these values , one has 12 equations in 14 unknowns : and V r of 

the forward and reverse shocks, and V, p, p, B of the plasma inside the 
interaction region. Hence there is no intrinsic reason why either the forward 
or the reverse shock should be stronger. The reason for the observed 
asymmetry must lie in the manner in which the shocks form. 

If one assumes that the formation of an infinite radial gradient in 
density is a sufficient condition for the formation of a shock, then one can 
attempt to estimate the heliocentric distance at which shocks will form. 

Figure 12a shows a stream seen by Voyager 1 between day 288 and day 300 
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of 1977 when Voyager 1 was at 1.2 AU. In this stream the density is 
aproximately anti-correlated with the bulk velocity. This rough 
anti-correlation of density and bulk velocity, in which the mass flux stays 
almost constant with respect to time, is typical of the flow near 1 AU. The 
velocity minimum and density maximum are located on day 291, and the velocity 

t 

maximum and density minimum are located on day 294. The region between days 
291 and 294 will later become the 'interaction region'. In this region, the 
solar wind plasma is being, compressed by the negative radial velocity gradient. 

Note that the velocity profile of this stream is not completely 
anti-symmetric about the midpoint of the interaction, on day 292. The 
velocity gradient in the high speed stream, after the passage of the 
interaction region, is less than the velocity gradient in the low speed 
stream. This structure is typical [Gosling et al., 1976 1 . As a consequence, 
the solar wind number density in the low speed stream will decrease with 
heliocentric distance more rapidly than will the number density of the high 
speed stream. 

The structure of a typical stream at 1,7 AU, still prior to the onset of 
stream-generated shocks , is shown in Figure 12b. The velocity profile retains 
roughly the same structure as the velocity profile at 1.2 AU, though the 
region of negative velocity gradient is narrower. The velocity minimum is on 
day 347 and the velocity maximum is on day 349. However, compression within 
the interaction region has radically changed the structure of the density 
profile. The density now has a maximum on day 348, near the center of the 
Interaction region. Also, as mentioned in the discussion of the stream at 1.2 
AU, the differing velocity gradients in those portions of the high and low 
speed streams which border the interaction region have led to different rates 
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of change of the dencity, so that the densities of the high and low speed 

k 

streams are now roughly equal. 

It has been assumed that the question of whether the forward or the 
reverse shock forms first will be determined by which radial density gradient 
steepens faster; the radial density gradient in the region in which the 
forward shock is to form, or the radial density gradient in the region in 
which the reverse shock is to form. But it is possible to calculate the rate 
of change of the radial density gradients which bound the interaction region. 
The gradients of thermal and magnetic pressures are still much less than the 
gradient of the dynamic pressure, so the only important MHD equation is the 
continuity equation. Considering only the radial derivatives, this becomes 
the one-dimensional continuity equation: 


&<«*) ■ -™ 2 f 


(26) 


The rate of change X of the radial gradient of NIT is then approximately: 


_ D a(NR 2 ) 

“ Dt 9R 


a ( nr 2 ) av / 2 v3 2 v 

3R BR / 5R 7 


(27) 


If one the stream in Figure 12B is typical of the streams at 1.8 AU, and 
that the stream pattern is being convected outwards at 400 km-sec”^, then the 
average plasma parameters for a compression region at 1.8 AU are shown in 
Table 3. 

At 1.8 AU the density gradient in the region that will become the forward 
shock is already steeper than the density gradient in the region that will 
become the reverse shock, and is steepening faster, which suggests that the 
forward shock will form first. Note that the forward shock is expected to 



22 


form first only because of the particular Initial conditions observed in this 
data. Obviously it 1 b possible to construct velocity and density profiles 
vhlch lead to an early formation of reverse shocks. A more exact 
determination of the location of the onset of forward and reverse shocks would 
require a detailed numerical model of solar wind stream evolution. Models of 
this sort have been constructed by Hundhausen (l973al and Pizzo (1982I, but 
are beyond the scope of this paper. 

Evolution and Structure 


Siscoe (1976 | predicted that on the average shock surfaces would conform 
to the Parker spiral. Figure 13 shows a plot of the azimuthal angle of the 
shock normals versus heliocentric distance. The value expected for shock 
surfaces which conform to the Parker spiral is shown in the curved trace. The 
apparent temporal variation in shock azimuthal angle is too large to allow a 
ready determination of the variation of the average azimuthal angle with 
heliocentric distance. However Packer spiral angle falls within the range of 
observed azimuthal angles. 

Figure 14 shows a plot of the north-south 3hock angles versus 
heliocentric distance. The average north-south angle remains fairly constant 
at zero. The variations about this mean are of the order of ten degrees. 

The large deviations of the shock azimuthal and north-south angles from 
the expected directions could be taken to mean that the shock surface is not 
smooth but contains large scale ripple3. Alternatively, it is possible that 
the boundaries of the co-rotating Interaction regions are not one continuous 
shock surface but consist of a large number of smaller intersecting shock 
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surfaces. To accurately reconstruct the large scale structure of co-rotating 
shocks, it will be necessary to include data from more than one spacecraft. 


Shock Energy Balance and Contribution to Global Heating 

* 

A detailed study of the energy balance for a shock requires knowledge of 
the electron temperatures before and after that shock. Electron temperatures 
for three shocks were provided by Sittler (private communication) . The 
contribution of the alpha particles to the total pressure was ignored. The 
shocks in question were the shocks of 1507 day 13 1979* 2056 day 17 1979* and 
1920 day 176 1979* The pre and post shock parameters for these shocks are 
listed in Table 4. The electron distribution functions have been divided in 
to a thermal 'core* distribution, and a non-thermal 'halo' distribution, and 
the parameters for these two distributions listed separatly. 

For the shocks of day 17 and day 176, the agreement between the measured 
pressure ratio and the predicted pressure ratio is quite good. The measured 
and the predicted pressure ratio (see below) do not agree quite as well for 
the shock of day 13, however the uncertainty in the measured pressure ratio 
for this event is quite large, primarily due to the uncertainty in unshocked 
electron density. 

From Equations 1 and 3, one can obtain an expression for the ratio of 
thermal energy density across a shock: 


p„ N, 1 B-. - 

r ■ t 1+r i( 1 - ip) + j (Mr 1 ) 2 )! 

P 1 w 2 P 1 a lt 


where ^ = P^/H^ 


( 28 ) 


p x = N^T-jAbJ/Bx) 

This relation is used to obtain the predicted pressure ratios shown in Table 4. 
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For the shocks of day XT and day 176, the agreement of the measured and 
the predicted ratio of pre and post shock thermal energy density is veil 
within the uncertainties, For the. shock of day 13, the agreement 1 b leas 
good. However, there is a large uncertainty in the observed shock parameters 

i 

for this shock. 

It vas not possible to obtain electron data for all of the shocks seen by 
Voyager 1. However, it is still possible to estimate the magnitude of the 
heating of the solar wind protons due to shocks; 

One can define an average proton heating rate (9E/at) per shock; 

a 

(3E/dt) * <AexU. >A (29) 

8 X 

where <Ae> is the average Jump in proton thermal energy density across a shock 
and A is the area of the shock surface. Running twelve shock averages of <Ae> 
calculated from the proton thermal speed and number density are plotted versus 
heliocentric distance in Figure 15* 

The average number density n_ of shocks is simply the inverse of the 
product of the area of a shock surface and the radial separation between the 
shocks or: 

n a = where <v> * average shock frequency (30 ) 

s 

The energy source term h g due to all shocks will then be the average 
heating rate per shock times the number density of the shocks, or; 

<v><Ae><U 1 > 

h s * — T v ~ > ~ ~ 
s 


( 31 ) 
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Typical values for <v>, <U.j>, and <V fl > for forward and reverse shocks, at 4 
and Q AU are shown in Table 5 , 

The energy source term due to shocks can be fit by a power law: 

<h Q R^> « h 0 R‘ P (where R is measured in AU) (38 ) 

h 0 »0.38*10" 14 erg-cnr^sec - ^, p*3«3 for forward shocks 

hQ*1.35 x 10"'^^ erg-cm“^sec‘' 1 l p»4.8 for reverse shocks. 

This energy source term is equivalent to the presence of a radial heat 
nux where: 

S S A 

- -<h/> 

Equation 33 can be solved to give: 

q (eff) R 2 = ^ R -a yhere ^ , .(! w ) J, and a = p _i 
-2 —2 -1 

q^«2. 47*10 erg-cm sec a=2.3 for forward shockB 

—2 —2 —1 

qQ=5*31 x 10 erg-cm sec a=»3*8 for reverse shocks. 

Equation 34 is an expression for the effective energy flux due to heating 
of the solar wind by stream-produced shocks. This expresssion is not valid at 
heliocentric distances less than 3 AU since stream-produced shocks do not 
appear closer to the sun. At 4 AU and 10 AU this effective energy flux is: 


(33) 


(34) 



2 6 


* 1* AU 
forward shockns q^ efl ^R 2 * l,02*lCf 3 
reverse shocks: q^ eff ^R 2 * 0.27*10" 3 


10 AU 



1.24xlO" 4 

erg cm" 2 

sec** 1 

0.84xl0" 5 

-2 

erg cm 

sec"' 1 ’ 


One can compare these 
the observed radial trends 


quantities with the heat 
in the solar wind energy 


flux necessary to explain 
fluxes, [Gazis, 1984): 


1 AU 4 AU IQ AU 

q R R 2 » 2.5«10“ 2 13.6xlO“ 3 6.9 x 10" 3 erg cm" 2 sec" 1 

At 4 and 10 AU, the effective energy flux due to heating of the 

solar wind by shocks is an order of magnitude less than the heat flux. The 

effective energy flux q^ eff ^ decreases slightly faster than q R and by 10 AU, 
(ef f ) 

q^ is a factor of 30 less than q R . 

One source of error in the determination of q^ e *^ is the failure to 

s 

include shocks with velocity jumps which lie below the detection threshold of 

20 km-sec . This omission reduces the measured value of '. The other 

major source of error is uncertainty in the determination of the power law 

exponent (3 for the variation of h with heliocentric distance. This 

8 

► / 

introduces an uncertainty in the magnitude of q^ e ' of at least a factor of 

8 

two. 

The net result of heating the solar wind protons by interplanetary shocks 

/ 

is the addition of an effective energy flux q„ ' . This effective energy 

s 

flux is comparable to or less than the heat flux q R required to explain the 

observed radial profiles of proton and electron temperature. As with q R , 

(eff) 

q^ tends to zero as R tends to infinity. 
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Quasi-Parallel Shocks 

A parallel shock should be characterized by a normal parallel to the 
magnetic field vector and by the fact that the average B-field magnitude does 
not change across the shock. 

There are seven shocks for which the angle between the shock normel and 

o 

the magnetic field vector is less than 30 and is less than the uncertainty in 
the normal uetermination , or for which the B-field magnitude does not vary by 
more than lOJt across the shock. Those shocks considered here to be 
quasi-parallel are listed in Table 6. There are additional shocks, most 
notably the shock of 02^7 day 296, 1978 which have structure similar to that 
which might be expected for parallel shocks and have sufficiently high 
uncertainties in some basic parameters that they might be qua si-parallel 
shocks. However the uncertainties make these Bhocks hard to analyze and they 
have been excluded from the table. 

The seven shocks listed in Table 6 have a distinctive structure: 

1) Each of these events has sign. 'leant wave activity upstream and 
downstream of the shock. In some cases, the wave amplitude 1 b as high as 30^ 
of the shock amplitude. This structure is similar to the structure usually 
seen upstream of Earth's bow shock in the dawn-to-noon quadrant. 

2) Each of these events has a finite thickness, of the order of ten 

5 

minutes. This corresponds to a spatial thickness of order of 3*10 km. 
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Slow Shock? 

As was diseased earlier, slow shocks should be characterised by M^ n <l and 
a Jump in B-field magnitude of opposite Blgn to the Jump in density. Eleven 

of the shocks observed at Voyager 1 fulfill these criteria, and they are 

» 

listed in Table 7» 

Both forward and reverse slow shocks are observed, though as for fast 

shocks, the forward shocks are the more numerous. All of these events are 

quite weak, with velocity Jumps less than 40 km-sec"' L and density ratios less 

than two. The orientations of these shocks are unusual; the shock normals are 

o 

directed in some cases as much as 90 away from the radial direction. 

Finally, the slow shocks show a tendency to occur very close together in 

time; six of the eleven events observed had some neighboring slow shock occur 

» 

within eight days. This suggests that some process may exist which generates 
strings of Blow shocks. 

Some of these shocks appear to be associated with shock pairs. If bo, 
and if the fast and slow shocks were formed simultaneously, their separation 
might provide a measure of how long ago these shocks were formed. 

Summary and Conclusions 

Between late 1977 and late 1980, the principle source of interplanetary 
shocks observed by the Voyager 1 spacecraft at large heliocentric distances 
was the interaction of high and low speed streams in the solar wind. These 
'stream-produced shocks' begin to appear at a heliocentric distance of between 
2 and 3 AU and persist out to the limits of these observations. 
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Forward shocks appear earlier than reverse shocks in these observations’, 
between 2 to 3 AU as opposed to 4 AU. The forward shocks are also 'stronger' 
than the reverse shocks; The forward shocks have a higher average speed with 
respect to the solar wind, a higher average velocity Jump, and a higher 
average density ratio than the reverse shocks. 

These stream-produced shocks appear to decay as they propagate outwards 
from the sun. This decay is. consistent with a simple model of shock evolution 
but it is not clear, from observations from only a single spacecraft, whether 
this decay represents a true radial variation, whether it is due to temporal 
effects or whether it is due to selection, effects in the shock survey. 

It is worth noting that stream-generated shocks can persist even after 
the streams that generated them have eroded away. Thus, even if the solar 
wind at large heliocentric distances is dominated by 'pressure waves' rather 
than streams as has been suggested by Burlaga (19831, one would still expect 
to observe stream-generated shocks. 

The majority (70^) of the quasi-perpendicular shocks are fast shocks. 

The Alfvenic mach numbers range from 1 to 9 and the density ratios range from 
1.5 to 4. There are also 11 cases of perpendicular slow shocks, and 7 cases 
of quasi-parallel shocks. 

For three shocks there are corresponding electron measurements. The 
measured and the predicted pressure ratios across the shocks agree, implying 
that there is no need to invoke the acceleration of high energy particles to 
explain the energy balance of these events. However the uncertainty in the 
pressure measurements is sufficiently large that as much as 10^ of the total 
energy density as measured in the shock frame could be lost into high energy 
particles or waves without being noticable. 
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The effect of stream-produced shocks in heating the solar wind at large 
heliocentric distances can be expressed by the Introduction of an additional 
heat Bource into the MHD equations. At 4 AU, the effect of this heat source 
is comparable to or less than the heat flux required to explain the observed 
radial variation of the solar wind temperature |Gazis, 1984]. The effect of 
shock heating appears to decrease more rapidly vith heliocentric distance than 
does the heat flux. However it is not clear how much of this decrease is due 
to real radial variation, how much is due to temporal variations in the solar 
wind stream structure, and what the effect is of shocks with velocity Jumps to 
small to be detected by this survey. 
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Contract 953733 (JPL). 



31 


References 

Abraham-Shrauner B. and 3. H. Yun, Interplanetary shocks seen by the Ames 
plasma probe on Pioneer 6 and 7» J* Geophys. Res., 8l , 2097 » 1976 

Bavassano, B. , P. Marian! » and N. F. Ness, Pioneer 8 Observations and 
interpretations sixteen interplanetary shock waves observed in 1968, J. 
Geophys. Res., 78 , 4535 » 1973 

Borrini, G. , J. T. Gosling, S. J. Bame, and W. C. Feldman, An analysis of 
shock wave disturbances observed at 1 AU from 1971 through 1978, J. Geophys. 
Res., 81 , 4363, 1982 

Burlaga, L. F. , and J. K. Chao, Reverse and forward slow shocks in the solar 
wind, J « Geophys , Res 76 , 7516, 1971 

Burlaga, L. F. Corotating pressure waves without fast streams in the solar 
wind, J. Geophys. Res., 88 , 6085, 1983 

Burlaga, L. F. , A., J. Lazarus, R. Schwenn, and H. Rosenbauer, Dynamical 
evolution of interplanetary magnetic fields and flows between 0.3 and 8.5 
AU, Geophys. Res. Lett. , 10 , 4l3, 1983. 

Chao, J. K. , and R. H. Lepping, A correlative study of SSC's, interplanetary 
shocks, and solar activity. J. Geophys. Res., 79 , 1799, 1974 

Colburn, D. S., and Sonnet, C. P. , Discontinuities in the solar wind. Space 
Scl. Rev., 5 . 439, 1966 

Dryer, M. , Z. K. Smith, G. Endrad, and J. H. Wolfe, Pioneer 7 observations of 
the August 29, 1966 interplanetary shock-wave ensemble, Cosmic 
Electrodynamics, 3 , 184, 1972 

Dryer, M. , Z. K. Smith, T. Unti, J. D. Mihalov, B. F. Smith, J. H. Wolfe, D. 

S. Colburn, and C, P. Sonett, Pioneer 9 and Ogo 5 observations of the 



32 


interplanetary multiple shock ensemble of February 2 1969, J. Geophys. Res., 
80, 3225, 1975 

Dryer, M., Z. K. Smith, E. J. Smith, J. D. Mlhalov, J. H. Wolfe, R. S. 

Steinolfsen, and S. T. Wu, Dynamic modeling of solar wind co-rotating stream 
interaction regions observed by Pioneer 10 and 11, J. Geophys. Res., 83 . 
4347, 1978 

Gazis, P. R. and A. J. Lazarus, The radial evolution of the solar wind, 1-10 
AU, Solar Wind 5 , 509, NASA Conference Publication 2280, ed, M. Neugebayer 
Gazis, P. R., Observations of plasma bulk parameters and the energy balance of 
the solar wind between 1 and 10 AU, J. Geophys, Res., 89 . 775, 1984 
Gosling, J. T., A. J. Hundhausen, and S. T. Bame, Solar wind stream evolution 
at large heliocentric distances: experimental demonstration and test of a 
model, J. Geophys, Res,, 81 . 2111, 1976 
Hundhausen, A. J,, and R. A. Gentry, Numerical solutions of flare-generated 
disturbances in the solar wind, J. Geophys. Res., 74 , 2908, 1969 
Hundhausen, A. J,, S. J. Bame, and M. D. Montgomery, Large scale 

characteristics of flare associated solar wind disturbances, J. Geophys. 
Res., 75 , 4631, 1970 

.Hundhausen, A. J., Nonlinear model of high-speed solar wind structure beyond 1 
AU, J. Geophys. Res., 78 , 1528, 1973a 

Hundhausen, A. J., Evolution of large-scale solar wind strutures beyond 1 AU, 
J. Geophys. Res., 78 , 2035, 1973b 

Lepping, R. P. , and Argentiero, P. D., Single spacecraft method of estimating 
shock normals, J. Geophys. Res., 76 , 4349, 1971 
Mlhalov, J. D. and J. H. Wolfe, Pioneer 10 studies of interplanetary shocks at 
large heliocentric distances, Geophys. Res. Lett., 6 , 491, 1979 



33 


Parker, E. N. Interplanetary Dynamical Processes , pl52, Intersclence, N. Y. 
1963 

Pizzo, V., A three dimensional model of corotating streams in the solar wind. 

3. Magnetohydrodynamic streams, J. Geophys. Rea., 87 . 4347, 1902 
Rosenau P. t and S. T. Suess, Slow shocks in the interplanetary medium, J. 
Geophys. Res., 82 , 3649. 1977 

Sarabai, V., Some consequences of non-uniformity of solar wind velocity, J . 
Geophys. Res., 68 , 1555 » 1903 

Siscoe, G. t Three dimensional aspects of interplanetary shock waves, J . 
Geophys. Res., 8l . 6235* 1970 

Smith, E. J., and J. H. Wolfe, Observations of interaction regions and 
co-rotating shocks between one and five AUs Pioneers 10 and 11, Geophys . 

Res. Lett.. 3 . 137, 1970 

Smith, E. J., and J. H. Wolfe, Observations of evolving solar wind streams and 
shocks beyond 1 AU, Study of Travelling Interplanetary Phenomena/ 1977 . Ed. 
M. A. Shea, et al., 227, D. Reidal, Hingham, Mass., 1977 
Smith, E. J., and Wolfe, J. H., Fields and plasmas in the outer solar system, 
Space Sci. Rev., 23 , 217, 1979 

Sonnet, C. P«, D. S. Colburn, L. Davis Jr., E. J. Smith, and P. J. Colman Jr., 
Evidence for a collision-free magnetohydrodynamic shock in interplanetary 
space, Physs. Rev. Letters. 13 . 153, 1904 
Steinolfsen, R. S.» M. Dryer, and Y. Nakagawa, Interplanetary shock pair 
disturbances: comparison of theory with space probe data, J ♦ Geophys . Res . , 
80, 1989, 1975 

Steinolfsen, R. S., and M. Dryer, Numerical simulations of MHD shock waves in 
the solar wind, J. Geophys. Res., 83 . 1570, 1978 



34 


Figure 1: Schematic representation of the structure of a co-rotatlng 

interaction region and its accompanying shocks, as viewed from above the plane 
of the solar equator. 

Figure fa Spacecraft trajectories for Voyager 1 and 2. 

Figure 3: Histogram of the shock velocity Jumps. 

Figure 4: Time axis plot of various solar wind bulk parameters observed 

during the shock event of 0455 day 313 1979* 

Figure 5* Histogram of the angle between the shock normals calculated 
from the mass continuity equation alone and the shock normal calculated 
including the effect of the B- field for the 71 fast shocks. 

Figure 6: Histogram of the measured density ratio divided by the density 

ratio predicted from the tangential pressure equation for all 83 shocks. 

Figure 7: Twelve shock running averages of shock frequency plotted 

versus heliocentric distance for forward shocks and reverse shocks. There is 
a one shock slip between the successive averages. 

Figure 8: Twelve shock running averages of shock speed plotted versus 

heliocentric distance for forward shocks and reverse shocks. There is a one 
shock slip between the successive averages. 

Figure 9: Twelve shock running averages of shock velocity Jump plotted 

versus heliocentric distance for forward shocks and reverse shocks. There is 
a one Bhock slip between the successive averages. 

Figure 10: Twelve shock running averages of shock density ratio plotted 

versus heliocentric distance for forward shocks and reverse shocks. There is 
a one shock slip between the successive averages. 

Figure 11: Schematic representation of plasma parameters in the shock 
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frame and in the spacecraft frame for a forward shock and for a reverse shock. 

Figure 12: Time axis plot of solar wind bulk velocity and density for a 
stream observed at 1.2 AU and a stream observed at 1.8 AU. The interaction 
regions are marked. 

Figure 13: The individual values of azimuthal angle of the Bhock normals 

for all 83 analyzed shocks plotted versus heliocentric distance. The solid 
trace shows the expected value if the shock surfaces were to conform exactly 
to the Parker spiral. 

Figure lU : The individual values of north-south angle of the shock 

normals for all 83 analyzed shocks plotted versus heliocentric distance. 

Figure 15: Twelve shock running averages of the Jump in energy density 

across each shock plotted versus heliocentric distance all shocks. There is a 
one shock slip between the successive averages. 
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Table 1. Voyager 1 Shock Times 


Year 

Day 

Hour 

Type 

Year Day Hour 

Type 

Year Day Hour 

Type 

1977 

301 

0655 

F 

1978 

( 287-288 

F) 

1979 42 0557 

F 


331 

2221 

F 


289 

2205 

R 

56 0308 

R 


347 

1215 

F 

■ 

296 

.0247 

G 

1 

Jupiter encounter 59-82 

(1978 

33 

1650 

F 


300 

1553-1620 R 

92 1139 

F 


178 

2303 

F 


305 

1314 

G 

93 0716 

R 

204 2258 

F 



316 1839-2255 F 

112 1140-1924 

R 


209 

1312 

F 


321 

0828 

S 

117 1712 

R 


213 

l44o 

F 


338 

0122 

G 

118 0430 

F 


217 

0220 

S 


338 1016-1324 F 

128 0520 

F 


223 

0825 

F 


341 0946-1343 R 

133 0244 

F 


228 0614-1537 R 


(356-357 

F) 

133 1057-1812 R 


237 1349-1942 F 


360 

0432 

F 

139 2045 

F 


24l 1543-1944 R 

1979 

1 

0628 

G 

152 0229 

F 


(251 

1938 

F) 


13 

1507 

F 

l 6 l 0235 

G 


264 0229-1006 F 


17 

2057 

R 

( 166-167 

F) 


267 

2142 

R 


21 

1334 

G 

(174-175 

F) 


279. 

1334 

R 


29 

1850 

G 

176 1920 

F 


284 

2130 

R 


37 

2206 

F 

179 0205 

R 


F = forward shock 
G = forward slope 
R = reverse shock 
S = reverse slope 

unconfirmed events are marked with parenthesis 
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Table 

1» Voyager 1 Shock Times (cont) 



Year Day 

Hour 

Type 

Year Day Hour 

Type 

Year Day Hour 

Type 

1979 191 

0946 

F 

1980 . 14 0736 

R 

1980 177 

0777 

F 

197 

1850 

F . 

28 0746 

F 

178 

0631 

F 

(206 

1630 

F) 

32 0616 

R 

201 

1511 

F 

211 1817-2150 R 

34 2246 

R 

222 

1955 

F 

217 

1120 

F 

(51 0500 

G) 

260 

2112 

S 

220 

1714-2132 R 

57 0122 

R 

261 

2010 

G 

222 0641-1539 R 

(60-61 

F) 

262 

0900 

G 

226 

0531 

R 

(82-83 

F) 




260-261 

R flare 

103 0305-0837 R 




266 

0308 

R 

112 1118 

F 




308 

2100 

G 

117 1447-2008 F 




313 

0455 

F 

122 1118 

F 




334 

0645 

G 

124 0828-1339 F 




339 

0134 

F 

146 0439 

F 




346 

0403 

F 

146 0817 

F 




346 

1800 

G 

150 0312 

F 




I960 7 

0936 

F 

(152 1320 

R) 




11 0536-1216 R 

(159 1203-1900 R) 

. 




F » forward shook 
G = forward slope 
R = reverse shock 
S a* reverse slope 

Unconfirmed events are marked with parenthesis 
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l'ABLE 2. Typical Values and Rates of Change for Shock Parameters 
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TABLE 3. Typical Solar Wind Parameters at 1.8 AU, Prior to Shock Formation 


Forward Side 


Reverse Side 


« 1 AU 2 


1.8*10 2 cm“ 3 AU -1 


-0.9X10 2 cm" 3 AU“ X 


4 v > 


a 6l0 km-8ec“ 1 'AU" X -6o6 km-sec’* X AU" X 


<-^5> 

1 AU 41 
*2 


?V> 


5R 

X 


1 au £ 


10 cm' 


-3 


10 cm 


-3 


2.6*10 3 km-ssec" 1 AU” 2 -2,6xi0 3 km-sec” 1 AU“ 2 


5.4*10“ 16 < an" 1 *8ec“ 1 


ts 


3.2xi0 -1 ^ cnf^sec" 1 




TABLE 4- Energy Balance for Three 8hocks 


Uo 


Time 1979 13 1507 1979 IT 2056 1979 176 1920 


upstream downstream 

upstream downstream 

upstream downstream 

k p 

0.073 

0.250 

0.150 

0.240 

0.715 

1.500 

N c 0.0581.012 0.2501.020 

0.1501.010 

0.2501.020 

1.6001.100 

0.9701.010 

N h 0.0l8t.004 0 

.0231.004 

0.0151.003 

0.0201.007 

0.0101.003 

0.0081.001 

t^io 1 *) 

O.67 

7.4 

2.2 

4.1 

3.5 

5.0 

T C (10 4 ) 

1.7±0.2 

6.810.2 

3.110.2 

4.010.4 

3.310.1 

3.510.2 

T h (l0 4 ) 

3211 

55120 

2615 

4313 

63113 

58117 

c 

iH 

D 

152 

49 

56 

33.6 

38 

18 

B lt 

0.71 

1.66 

0.34 

0.60 

0.58 

O.89 

p 1 ><10 13 

0.67±.07 

2.61.03 

0.461.03 

1.4+.10 

3.5±.70 

10 .'4+ .30 

P c *10 13 

1.41.30 

2.31.10 

o.64i.o4 

i.4i.io 

4.41.10 

7.71.50 

Ph^lO 13 

0. 81.20 

1.81.80 

0.541.18 

1.21.30 

0.691.15 

0.801.20 

r l 

9.81.90 


4.81.50 


0.801.20 


1/Pi 

0.71.07 


0.281.02 


0.16 


(“) 

'■p^ 'meas 

2.3H.0 


2.410.2 


2.210.4 


'>p 1 ' l pred 

4. 8+0. 4 


2.310.2 


1.910.2 





TABLE 5. Shock Frequencies, Speeds, and Heating Rates 


hi 




Forward Shocks 

Reverse ShockB 



4 AU 

8 AU 

4 AU 

8 AU 

<v> 

S3 

I.l6xl0 -6 

0.58xl0" 6 

g 

0 . 69 x 10 

» 

6 1 
0.23x10 sec 

«V 

& 

160 

0 

CO 

120 

-1 

50 km sec 

<v 

X 

570 

490 

290 

360 km sec" 

< 4 '<rsj )2 > 

3 

1.2xi<f 10 

4.2X10“ 9 

6.2X10" 9 

2.0X10" 9 erg-cnf 3 


a 

3.9*10" 17 

h.OxicT 18 

l.TTxiO*' 17 

6.5X10”’ 3,9 erg-cm^sec” 1 
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TABLE 6. Quasi-Parallel Shocks 


T 


Year 

Day 

Hour 

Angle 

<w 

1978 

20k 

2258 

0 

150 

1.0 


305 

1314 

' 152° 

1.5 


338 

0122 

140° 

0.9 

19T9 

5 6 

0508 

168 ° 

1.1 


260 

2329 a 

163° 

1.1 


2 66 

0308 

150° 

0.9 

1980 

261 

2010 

6 o° 

1.0 


a The time of his event 

la listed in 

Table 1 as 1979 260 - 261 . 


1*3 


TABLE 7. Slow Shocks 


Year Day Hour 

1979 21 1334 

29 1846 

92 1139 

93 0716 
197 1850 
334 0645 
346 0403 
346 1800 

1980 32 0616 
122 1118 
146 0439 
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Schematic Representation of Co-Rotating 
Interaction Regions 


Figure I 
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Reverse Shock 



Figure II 
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